Line Shapes of the X(3872) 



Eric Braaten and Meng Lu 

Physics Department, Ohio State University, Columbus, Ohio 43210, USA 

(Dated: February 1, 2008) 

Abstract 

If the quantum numbers of the X(3872) are J'^*^ = 1^^, the measurement of its mass implies 
that it is either a loosely-bound hadronic molecule whose constituents are a superposition of the 
charm mesons pairs D*^D^ and D^D*^ or else it is a virtual state of these charm mesons. Its 
binding energy is small enough that the decay width of a constituent D*^ or D*^ has a significant 
effect on the line shapes of the X resonance. We develop a simple approximation to the line 
shapes that takes into account the effect of the D*^ width as well as inelastic scattering channels 
of the charm mesons. We carry out a simultaneous fit to the line shapes in the J/'ijjTT^TT^ and 
D^D^TT^ channels measured in the decays — > + X hy the Belle Collaboration. The best 
fit corresponds to the X(3872) being a bound state just below the D*^D^ threshold, but a virtual 
state just above the D*^D^ threshold is not excluded. 

PACS numbers: 12.38.-t, 12.39. St, 13.20.Gd, 14.40. Gx 
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I. INTRODUCTION 



The X(3872) is a hadronic resonance near 3872 MeV discovered in 2003 by the Belle 
Collaboration [1] and subsequently confirmed by the CDF, Babar, and DO Collaborations 
[2-4]. In addition to the discovery decay mode J / t:^ Ti^ , the X has been observed to 
decay into J/ip'^ and J/V' vt+tt^tt^ [5]. The decay into J/ip^ implies that the X is even 
under charge conjugation. An analysis by the Belle Collaboration of the decays of X into 
J/ipTi~^Ti~ strongly favors the quantum numbers J^*^ = l"*""*", but does not exclude 2"'""'" [6]. 
An analysis by the CDF Collaboration of the decays of X into J/V'tt+tt" is compatible 
with the Belle constraints [7]. The Belle Collaboration has also discovered a near-threshold 
enhancement in the D^D^n^ system near 3875 MeV [8]. If this enhancement is associated 
with the X(3872), the tiny phase space available would rule out J = 2, leaving 1++ as the 
only option. 

An important feature of the X(3872) is that its mass is extremely close to the D*^D° 
threshold. The PDG value for Mx comes from combining measurements of X in the 
J/ipT:~^T:~ decay mode [9]. After taking into account a recent precision measurement of 
the D° mass by the CLEO Collaboration [10], the difference between the PDG value for Mx 
and the D*^D^ threshold is 



where M*o and Mq are the masses of D and D . The negative central value is compatible 
with the X being a bound state of the charm mesons. The PDG value for the mass of 
X(3872) comes from combining measurements of X in the J/ijjTr~^7r~ decay mode [9]. The 
peak of the near-threshold enhancement in £)°L)°7r° [8] is at a mass satisfying 



We have obtained this result from the value of M quoted in Ref. [8] by subtracting 2Mo -|- 
(M*o — -^o); where Mq is the 2006 PDG fitted mass for the D'^, and by dropping the error 
bar associated with the mass. The positive central value in Eq. (2) is compatible with 
X being a virtual state of charm mesons. The difference between the masses in Eqs. (2) and 
(1) is 4.7tJ 8 MeV. This is more than two standard deviations, which raises the question of 
whether the decays into J/ijj7r~^7i~ and D^D^tt^ are coming from the same resonance. 

The proximity of the mass of the X(3872) to the threshold has motivated its 

identification as a weakly-bound molccTilc whose constituents arc a superposition of the 
charm meson pairs Z)*'^Z)'^ and D^^D*^^ [11-14]. The establishment of the quantum numbers of 
the X(3872) as would make this conclusion almost unavoidable. The reason is that these 
quantum numbers allow S-wave couplings of the X to D*^D^ and D^D*^. Nonrelativistic 
quantum mechanics implies that a resonance in an S-wave channel near a 2-particle threshold 
has special universal features [15]. Because of the small energy gap between the resonance 
and the 2-particle threshold, there is a strong coupling between the resonance and the 
two particles. This strong coupling generates dynamically a large length scale that can be 
identified with the absolute value of the S-wave scattering length a of the two particles. 
Independent of the original mechanism for the resonance, the strong coupling transforms 
the resonance into a bound state just below the two-particle threshold if a > or into a 
virtual state just above the two-particle threshold if a < 0. If a > 0, the bound state has a 
molecular structure, with the particles having a large mean separation of order a. 



Mx - {M^Q + Mo) = -0.6 ± 0.6 MeV, 



(1) 



M - (M,o + Mo) = +4.1 ± Q.7tii MeV. 



(2) 
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The universality of few-body systems with a large scattering length has many applications 
in atomic, nuclear, and particle physics [16]. To see that these universal features are relevant 
to the X(3872), we need only note that its binding energy is small compared to the natural 
energy scale associated with pion exchange: m^/(2M*oo) ~ 10 MeV, where M*oo is the 
reduced mass of the two constituents. The universal features of the X(3872) were first 
exploited by Voloshin to describe its decays into D^D^n^ and D^D^'j, which can proceed 
through decay of the constituent D*^ or D*^ [14] . Universality has also been applied to the 
production process B — > KX [17, 18], to the hne shapes of the X [19], and to decays of X 
into J/t/j and pions [20]. These apphcations rely on factorization formulas that separate the 
length scale a from all the shorter distance scales of QCD [19]. The factorization formulas 
can be derived using the operator product expansion for a low-energy effective field theory 
[21]. 

Other interpretations of the X{3872) besides a charm meson molecule or a charm meson 
virtual state have been proposed, including a P-wave charmonium state or a tetraquark state. 
(For a review, see Ref. [22].) The discrepancy between the masses Mx and M in Eqs. (1) 
and (2) has been interpreted as evidence that the J/ ip tt+tt" events and the D^D^ti^ events 
arise from decays of two distinct tetraquark states whose masses differ by about 5 MeV 
[23] . If the charmonium or tetraquark models were extended to include the coupling of the 
X to D*^D^ and D^D*^ scattering states, quantum mechanics implies that the tuning of 
the binding energy to the threshold region would transform the state into a charm meson 
molecule or a virtual state of charm mesons. Any model of the X(3872) that does not 
take into account its strong coupling to charm meson scattering states should not be taken 
seriously. 

We will assume in the remainder of this paper that the quantum numbers of the X(3872) 
are 1"*"^, so that it has an S-wave coupling to D*^D^ and D^D*^. In this case, the measured 
mass AIx implies unambiguously that X must be either a charm meson molecule or a 
virtual state of charm mesons. The remaining challenge is to discriminate between these 
two possibilities. The Belle Collaboration has set an upper bound on the width of the X: 
Tx < 2.3 MeV at the 90% confidence level [1]. If the width of the X is sufficiently small, 
there is a clear qualitative difference in the line shapes of X between these two possibilities. 
We first consider the D^D'^n'^ decay mode, which has a contribution from the decay of a 
constituent D*'^. If the X was a charm meson molecule, its line shape in D^D^tt^ would 
consist of a Breit-Wigner resonance below the D*^D^ threshold and a threshold enhancement 
above the D*^D^ threshold. If the X was a virtual state, there would only be the threshold 
enhancement above the D*^D^ threshold. We next consider decay modes that have no 
contributions from the decay of a constituent such as J/^ tt+tt^. If the X was a charm 
meson molecule, its line shape in such a decay mode would be a Breit-Wigner resonance 
below the D*^D^ threshold. If the X was a virtual state, there would only be a cusp at 
the D*^D^ threshold. The possibility of interpreting the V(3872) as a cusp at the D*^D^ 
threshold has been suggested by Bugg [24] . Increasing the width of the X provides additional 
smearing of the line shapes. This makes the qualitative difference between the line shapes 
of a charm meson molecule and a virtual state less dramatic. To discriminate between these 
two possibilities therefore requires a quantitative analysis. 

The proposal that X(3872) is a charm meson virtual state has received support from 
a recent analysis by Hanhart et al. [25] of data on + J/i(j7r~^7r~ and B~^ — >■ 

+ D^D^Ti^ from the Belle and Babar Collaborations. They concluded that the D^D^tt^ 
threshold enhancement is compatible with the X(3872) only if it is a virtual state. One flaw 
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in the analysis of Ref. [25] is that they did not take into account the width of the constituent 
They identified the rate for D^D^n^ with the rate for D*^D^ and D^D*^ multiphed by 
the 62% branching fraction for D*^ — > D^tt^. Thus they assumed that D^D^tt^ could not be 
produced below the D*^D^ threshold. Since a bound state has energy below the threshold, 
this implies that X cannot decay into D^D^n^ if it is a bound state. This contradicts 
one of the universal features of an S-wave threshold resonance. As the binding energy of 
the resonance decreases, the mean separation of the constituents grows increasingly large, 
suppressing all decay modes except those from the decay of a constituent. 

In this paper, we develop a simple approximation to the line shapes of the X(3872) that 
takes into account the D*^ width. In Sec. II, we analyze the decays of the D* mesons and 
give a simple expression for the energy-dependent width of a virtual D*^ . In Sec. Ill, we 
give an expression for the scattering amplitude for D*^D^ that takes into account the width 
of the D*^ as well as inelastic scattering channels for the charm mesons. In Sec. IV, we use 
that scattering amplitude together with the factorization methods of Rcfs. [19, 21] to derive 
the line shapes of X(3872) in the decays + X. In Sec. V, we give an order-of- 

magnitude estimate of the short-distance factor associated with the transition in 

the factorization formula. In Sec. VI, we discuss the analysis of Ref. [25] and point out the 
flaw from neglecting the effects of the D*^ width. In Sec. VII, we carry out a new analysis of 
the data on B+ K+ + Jj^ tt+tt" and 5+ ^ 7^+ + D^D^ti^ from the Befle Collaboration, 
taking into account the effects of the D*^ width. In Sec. VIII, we discuss ways in which our 
description for the line shapes could be further improved. 



II. D* WIDTHS 

In this section, we present a quantitative analysis of the decay widths of the D* mesons. 
There are particles with six different masses that enter into the analysis. We therefore 
introduce concise notation for the masses of the charm mesons and the pions. We denote 
the masses of the spin-0 charm mesons D'^ and by Mq and Mi, respectively. We denote 
the masses of the spin-1 charm mesons D*^ and D*'^ by M*o and M*i, respectively. We 
denote the masses of the pions vr'' and tt"'" by and mi, respectively. In each case, the 
numerical subscript is the absolute value of the electric charge of the meson. The pion 
mass scale corresponding to either mo or mi will be denoted by rriT^. The result of a recent 
precision measurement of the mass by the CLEO Collaboration is Mq = 1864.85 ± 0.18 
MeV, where we have combined the errors in quadrature [10]. We use the PDG values for 
the differences between the charm meson masses [9]. The errors on the pion masses are 
negligible compared to those on the charm meson masses. It is also convenient to introduce 
concise notations for simple combinations of the masses. We denote the reduced mass of a 
spin-1 charm meson and a spin-0 charm meson by 

M., = (3) 
We denote the reduced mass of a pion and a spin-0 charm meson by 

m, + Mj- 

The phase space available for the decay D* — > Dn depends sensitively on the difference 



4 



between the D* mass and the Dn thresholds: 



5ijk = M^i - Mj - ruk. (5) 

The differences between the D* masses and the thresholds for Dti states with the same 
electric charge are 

(5ooo = 7.14 ± 0.07 MeV, (6a) 
(5oii = -2.23 ±0.12 MeV, (6b) 
5no = 5.66 ±0.10 McV, (6c) 
= 5.85 ± 0.01 MeV. (6d) 

The T-matrix elements for the decays D* Dn are proportional to the dot product of 
the polarization vector of the D* and the 3-momentum of the vr. The coefficient of the dot 
product can be expressed as the product of a constant, which we denote by ^^/3/2 g/ f^r, and 
a Clebsch-Gordan coefficient. The partial widths for the decays D* — > Dn are 

rp.+ ^ = '-ffm%'S%l (7a) 

rp-+ D-A = ^yXl (7b) 

rp-» ^ DV] = ^rnl'S^. (7c) 

Since ^on < 0, D*^ does not decay into D~^tt~. 

The PDG value for the total width of the D*+ is T[D*+] = 96±22 keV [9]. We can use the 
PDG values for the branching fractions for D*+ D+n^, D*+ D\+, and D*+ D+-f 
to determine the partial widths for those decays: 

r[D*+ D\+] = 65.0 ± 14.9 keV, (8a) 
r[D*+ D+7r°] = 29.5 ± 6.8 keV, (8b) 
r[D*+ D+7 ] = 1.5 ± 0.5 keV. (8c) 

By fitting the partial widths in Eqs. (8a) and (8b) to the expressions in Eqs. (7a) and (7b), 
we obtain consistent determinations of the constant g/ f.^^. The more accurate of the two 
determinations comes from the decay D*+ L'°7r+: g/U = (2.82 ± 0.32) x 10"^ MeV"^/^. 
Inserting this value into Eq. (7b), we can predict the partial width for D*^ D°7r°: 

r[D*o ^ D\°] = 40.5 ± 9.3 keV. (9) 

By combining this prediction with the PDG value for the branching fraction for — > D^tt^, 
we obtain a prediction for the total width of the D*°: r[D*^] = 65.5 ± 15.4 keV. We can then 
predict the partial width for the decay D*° — > D°7 using the PDG value for its branching 
fraction: 

r[D*o ^ D%] = 25.0 ± 6.2 keV. (10) 

Since the decay rates for D* — > Dtt in Eqs. (7) scale like the 3/2 power of the energy 
difference between the D* mass and the Dn threshold, they are fairly sensitive to the mass 
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FIG. 1: The energy-dependent width r*o(-E') of a virtual D*^ as a function of the energy E relative 
to the D*^ mass. The point with error bars at = indicates the value and uncertainty of the 
predicted width of 

of the D*. A virtual with energy M*o + E can be considered as a D*^ whose rest energy 
differs from its physical mass by the energy E. The width of the virtual particle varies with 
E. We denote the energy- dependent width of the D*'^ by T^o{E). li \E\ is small compared to 
the pion mass, the energy-dependent width can be obtained simply by scaling the physical 
partial width for the decay D°7r°: 

r,o{E) = T[D*' ^ + T[D*' ^ D\'] [ [{Sooo + E)/6ooof' ^^(^ooo + E) 

+2 {mn/moof [{Son + E)/6ooof' 0{6on + E)] . (11) 

We ignore any energy dependence of the decay widths into D^'-f, because the photon energy 
and the phase space for the decays D*^ D^'-f do not vary dramatically in the D*^D^ 
threshold region. In Fig. 1, we plot the energy-dependent width of the D*^ as a function of the 
energy E relative to the D*^ mass. The three terms in Eq. (11) have obvious interpretations 
as energy-dependent partial widths for decays of D*^ into D^'-f, D^'tt^, and D^n^ . We define 
energy-dependent branching fractions Btoq^{E), Brooo(-^)5 and Broii(-E) by dividing these 
terms by r^o{E). For example, the energy-dependent branching fraction for D*^ — > D^tt^ is 

Brooo(^) = 7^?°''°^ [(^000 + E)/5ooof' 0{5ooo + E). (12) 

III. LOW-ENERGY D*^D^ SCATTERING 

In this section, we discuss the low-energy scattering of the charm mesons D*^ and D^. 
The existence of the X(3872) with quantum numbers 1++ implies that there is an S-wave 
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resonance near threshold in the channel with even charge conjugation: 



{D*D)l = ^ + . (13) 

The isospin splittings between the charm meson masses are approximately 4.8 MeV for 
D+ - and 3.3 MeV for D*+ - D*^. The energy splitting between the D*+D- and D*^D^ 
thresholds is the sum of these isospin splittings: 

V = 8.08 ± 0.12 MeV. (14) 

We will use approximations that are valid when the energy E relative to the D*^D^ threshold 
is small compared to u. Thus the relative momenta of the charm mesons is required to be 
small compared to (2M,,ooi^)^^^ ~ 125 MeV. This is numerically comparable to the pion mass 
scale: fti 135 MeV. 

The presence of the X(3872) resonance so close to the D*^D^ threshold with quantum 
numbers that allow S-wave couplings to D*'^D^ and D'^D*'^ indicates that it is necessary to 
treat the interaction that produces the resonance nonperturbatively in order to take into 
account the constraints of unitarity. The resonance is in the channel (D*Z))*}_. There may 
also be scattering in the channel that is odd under charge conjugation, but we will assume 
that it can be neglected compared to scattering in the resonant channel. Since the kinetic 
energy is so low, the scattering will be predominantly S-wave. 

We express the transition amplitude A{E) for the scattering of nonrelativistically nor- 
malized charm mesons in the channel (D*i))1_ in the form 

-'"*00 

where f{E) is the conventional nonrelativistic scattering amplitude. An expression for the 
low-energy scattering amplitude with an S-wave threshold resonance that is compatible with 
unitarity is 

where k,{E) = (— 2M*oo-E' — isy/'^ and E is the total energy of the charm mesons relative to 
the D*^D'^ threshold. If E is real, the variable k{E) is real and positive for £^ < and it 
is pure imaginary with a negative imaginary part for E > Q. The parameter 7 in Eq. (16) 
can be identified as the inverse scattering length. If 7 is complex, the imaginary part of the 
scattering amplitude in Eq. (16) satisfies 

Im/(£;) = |/(^)|2lm [7 -/.(£;)]. (17) 

The scattering amplitude f{E) in Eq. (16) satisfies the constraints of unitarity for a 
single-channel system exactly provided 7 is a real analytic function of E. For positive real 
values of the energy E, Eq. (17) is simply the optical theorem for this single-channel system: 

Im/(£;) = |/(£;)|V2M,oo^ {E>Q). (18) 

The left side is the imaginary part of the T-matrix element for elastic scattering in the 
{D*Dy]^ channel. The right side is the cross section for elastic scattering multiplied by 
[£'/(2M*oo)]^^^- If 7 > 0, the amplitude f{E) has a pole at a negative value of the energy 
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E, indicating the existence of a stable bound state. If 7 varies sufficiently slowly with E, 
we can approximate it by a constant. The pole is then at Epoic ~ — 7^/(2M*oo) and the 
binding energy is 7^/(2M*oo)- In addition to the contribution to the imaginary part of f{E) 
in Eq. (18), there is a delta-function contribution at = Ep^ie' 

lmf{E) ^ ^(5(£; + 77(2M,oo)) (E < 0,7 > 0). (19) 

If 7 < 0, the pole in the amplitude f{E) is not on the real E axis, but on the second sheet 
of the complex variable E. The standard terminology for such a pole is a virtual state. The 
imaginary part of the amplitude is nonzero only in the positive E region and is given by 
Eq. (18). 

We would like to identify the bound state in the case 7 > with the X(3872). However, 
since the X decays, it must have a nonzero width. Its decay modes include D^D'^n^ and 
D^D'^j, which receive contributions from decays of the constituent D*^ or D*^. There are 
also other decay modes, including J/'j/'tt+tt", J/'?/^ 7r+7r~7r°, and J/ip'y. All the decay modes 
of X are inelastic scattering channels for the charm mesons. The dominant effects of these 
inelastic scattering channels on scattering in the (D*Z))° channel can be taken into account 
through simple modifications of the variables 7 and n{E) in the resonant amplitude f{E) in 
Eq. (16). The effects of the decays of the constituent D*'^ or D*^ can be taken into account 
simply by replacing the mass M^,o that is implicit in the energy E measured from the D*^D^ 
threshold by M*o — ir*o(-E)/2, where T^,q{E) is the energy-dependent width of the D*^ given 
in Eq. (11). This changes the energy variable k{E) defined after Eq. (16) to 

n{E) = V-2M,oo[£; + ir.o(£^)/2]. (20) 

At the threshold E — the energy-dependent width V^q{E) reduces to the physical width 
r[£)*°]. The expression for k{E) in Eq. (20) requires a choice of branch cut for the square 
root. If E is real, an explicit expression for k.{E) that corresponds to the appropriate choice 
of branch cut can be obtained by using the identity 



^J-2M[E + iT/2] = \fM 



+ r2/4 - - i (^^E^ + rV4 + E^^ 



(21) 



The effects of inelastic channels other than D^D^n^ and D^D^^ on scattering in the {D*D)'^ 
channel can be taken into account by replacing the real parameter 7 in Eq. (16) by a complex 
parameter with a positive imaginary part. With these modifications of the variables 7 and 
k{E). the expression for the imaginary part of the amplitude f{E) in Eq. (17) can now 
be interpreted as the optical theorem for this multi-channel system. The right side can 
be interpreted as the total cross section for scattering in the {D*D)^ channel multiplied 
by [£'/(2M*oo)]^^^- The term proportional to Im7 can be interpreted as the contribution 
from the inelastic scattering channels. This interpretation requires Im7 > 0. The term 
proportional to —Im «;(£') includes a contribution that reduces to the right side of Eq. (18) 
in the limits Im7 — > and r*o — > 0. It can be interpreted as due to elastic scattering 
into D*^D^ or D^D*^. If Re 7 > 0, the term proportional to —Im «;(£') also includes a 
contribution that reduces to the right side of Eq. (19) in the limits Im7 — and F^o ~^ 0. 
It can interpreted as due to scattering into the X resonance. The D*^ or Z)*° produced by 
the elastic scattering process will eventually decay. Similarly, the constituent D*^ or D*^ in 
the X resonance will eventually decay. Thus the ultimate final states corresponding to the 
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term proportional to —lmK,{E) in Eq. (17) can be identified as D^D^ti^ and D^D^^, and 
also D^D^Ti~ and D^D^n'^ if the energy E exceeds the threshold |5oii| = 2.2 MeV. 

The scattering amplitude f{E) given by Eqs. (16) and (20) is a double- valued function of 
the complex energy E with a square-root branch point and a pole. If we neglect the energy 
dependence of the width r*o(-E'), the branch point is near E — —ir[D*^]/2 and the position 
of the pole is 

If Re 7 > 0, the pole is on the physical sheet of the energy E. It can be expressed in the 
form 

£;poie = -Ex - irx/2, (23) 

where Ex and Fx are given by 

Ex ~ [(Re7)'-(Im7)2]/(2M,oo), (24a) 
Fx ~ r[L'*°]+2(Re7)(Im7)/M,oo. (24b) 

If Ex > and Fx -C Ex, the state X is a resonance whose line shape in the region 
\E — Ex\ <^ Ex is a Breit-Wigner resonance centered at energy —Ex with full width at half 
maximum Fx- If Fx is not small compared to Ex, one can choose to define the binding 
energy and the width of X to be given by the expressions for Ex and Fx in Eqs. (24), but 
they should not be interpreted literally. If Re 7 < 0, the pole at the energy Epoie in Eq. (22) 
is on the second sheet of the energy E and it corresponds to a virtual state. In this case, 
the expressions for Ex and Fx in Eqs. (24) have no simple physical interpretations. 



IV. LINE SHAPES OF X(3872) IN B+ DECAY 



If a set of particles C has total quantum numbers that are compatible with those of the 
X(3872) resonance and if the total energy E of these particles is near the D*^D^ threshold, 
then there will be a resonance in the channel C. The line shape of X(3872) in the channel 
C is the differential rate as a function of the total energy E of the particles in C . In this 
section, we discuss the line shapes of the X(3872) for energy E close enough to the D*^D^ 
threshold that we need only consider the resonant channel {D*D)^ defined in Eq. (13). This 
requires to be small compared to the energy u 8.1 MeV of the D*^D~ threshold. 

In Ref. [19], it was pointed out that the hne shape for decay modes of X that do not 
involve decays of a constituent can be factored into short-distance factors that are insensitive 
to E and to the inverse scattering length 7 and a resonance factor that depends dramatically 
on E and 7. In Ref. [21], it was shown that the factorization formulas could be derived using 
the operator product expansion for an effective field theory that describes the cc sector of 
QCD near the D*^D^ threshold. The factorization formulas hold for any short- distance 
process, which we define to be one in which all the particles in the initial state and all the 
particles in the final state beside the resonating particles have momenta in the resonance 
rest frame that are of order or larger. An example of a short-distance process is the 
discovery mode S+ — > K'^ + J/t/jn'^n'. In the J/t/) tt+tt" rest frame, the momenta of the 5+ 
and arc 1555 MeV, which is much larger than m,^. The root-mean-square momentum 
of the J/ip in the J/ipn'^ir' rest frame is 319 MeV, which is significantly larger than m^r. 

We consider the short- distance process — > -|- C in which the collection of particles 
denoted by C have a resonant enhancement associated with the X(3872). If the total energy 
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E of the particles in C is near the D*^D^ threshold, the amplitude for the process 

+ C factors into a short-distance factor associated with the process — > + {D* D) 
a resonance factor f{E) given by Eq. (16), and a short-distance factor associated with 
the process {D*D)^_^_ C. Since the short-distance factors are insensitive to E, the only 
dramatic dependence on E comes from the factor f{E). Thus the line shape dV/dE is 
proportional to \f{E)\'^. 

For an alternative derivation of this line shape that also gives the line shapes for D^D^tt^ 
and D^D^j, we consider the inclusive differential decay rate into all channels that are en- 
hanced by the {D*D)^ resonance. Using cutting rules, the inclusive decay rate in the region 
near the D*^D^ threshold summed over all resonant channels can be expressed in the form 

— [B+ -^K+ + resonant] = 2 Im / (£;) , (25) 
dE 

where is a short-distance factor defined in Ref. [21]. The optical theorem in Eq. (17) 
can be used to resolve this inclusive resonant rate into two terms. We interpret the term 
proportional to Im 7 as the sum of all contributions from short- distance channels. Thus the 
fine shape of X(3872) in a specific short-distance channel C can be expressed as 

^K+ + c] = 2rf: |/(E)|2r^(E). (26) 

where is the same short-distance factor as in Eq. (25) and V^{E) is a short-distance 
factor associated with the transition of the charm mesons to the particles in C. The factor 
r*^(i?) in Eq. (26) differs by a factor of vr from the factor F*^ defined in Ref. [21]. With 
this convention, T'~'{E) is the contribution of the state C to the term Im7 in Eq. (17). As 
indicated by the argument we have allowed for the possibility that the dependence of this 
term on E is not negligible in the D*^D^ threshold region. The justification for Eq. (26) 
relies on first order perturbation theory in the coupling to the channel C. If that coupling is 
too large, the channels {D*D)^ and C would have to be treated as a two-channel resonating 
system. In the case of the channel J/ipTT~^7T~ , first order perturbation theory can be justified 
by the small branching ratio of the decay of X into J/-^ tt+tt" relative to D^D^n'^ measured 
by the Belle Collaboration [8] . 

We interpret the term proportional to — Im k{E) in Eq. (17) as the sum of all contribu- 
tions from channels that correspond to D*^D^ or D^D*^ followed by the decay of the 
or D*^. We can further resolve this term into the contributions from the channels D^D^ir'^, 
D^D'^j, D'^D^TT~, and D^D~tt'^ by multiplying it by the energy-dependent branching frac- 
tions Brooo(-E'), Broo7(-E)! |Broii(i?), and |Broii(£'), which add up to 1. A simple expression 
for —lmK,{E) can be obtained by using the identity in Eq. (21). The resulting expression 
for the line shape of X in the D^D^tt^ channel is 

_ , r , 1 1/2 

_[5+ ^ K+ + D'D\'] = 2Ff: \f{E)\' [M,oo{VE' + r,o{E)yA + E)\ Brooo(^), 

(27) 

where F^+ is the same short-distance factor as in Eq. (25) and Brooo(-E') is given in Eq. (12). 

In Ref. [20], the decay rates of X into J/i/j plus pions and photons were calculated under 
the assumption that these decays proceed through couplings of the X to J/ip and the vector 
mesons p° and ui. The results of Ref. [20] can be used to calculate the dependence of the factor 
r^{E) in Eq. (26) on the energy E for C ^ J/V'TT+Tr", J/V' 7r+7r-7r°, J/V'7r°7, and J/V'7. 
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FIG. 2: Energy dependence of the final-state factors T'-''{E) for the channels C = J/iIj'k~^'k~ 
and J / Tp TT^ 7:~ TT^ . The factors r'-"(£')/r'-^(0) as functions of the energy E relative to the D*^D^ 
threshold are shown as solid curves. The dashed curves are the corresponding factors used in 
Ref. [25]. 



The normalization of this factor, which we can take to be r'^(O), can only be determined by 
measurements of the branching fraction of X(3872) into the final state C. In Fig. 2, we plot 
the energy dependence of the factors T'-^{E) for C = J/ip7i^TT~ and J/ip 7t~^t:~tt'^. The larger 
variations for J/ipTr~^7r~7r^ are due to the width of cu being much smaller than that of p°. 
Very close to the threshold, the final state factors can be approximated by the expressions 

pj/V-TT+TT-^^) ^ (l + a2^)Ff", (28a) 
pj/V-TT+Tr-TrO^^) ^ (l + a3^)rf", (28b) 

where the coefficients are a2 — 0.0175 McV"^ and ag = 0.0809 MeV"\ The approximations 
in Eqs. (28a) and (28b) are accurate to within 1% for -8.5 MeV < E < 10.1 MeV and 
-4.3 MeV < E < 1.0 MeV, respectively. 

In Fig. 2, we also show the energy dependence of the factors analogous to T^{E) that 
were used in Ref. [25]. Those factors were denoted by r^+^r-j/^ (-£)/(? and r^+^-^oj/^(£')/g'. 
In Ref. [25], the 7i~^7r~ and tt+tt"??'' systems were treated simply as Breit-Wigner resonances 
with the masses and widths of the and uj. The additional dependence on E from the 
coupling of the resonances to pions and from the integration over the phase space of the 
pious was not take into account. In the region jE'l < 1 MeV, the factors in Ref. [25] analogous 
to r'^"(i!^)/r*^(0) differ from ours by less than 1.4% for C = J/ijj 7r~^'K~ and by less than 2.4% 
for C = J / ip 71^ 7r~ 71^ . As illustrated in Fig. 2, the differences are more substantial when \E\ 
is 10 MeV or larger and they are particularly large for J/tjj Tr'^ir^iT^ . 

In Fig. 3, we illustrate the hne shapes of X (3872) in the D'^D^n^ channel and in a short- 
distance channel, such as J/'07r'^7r~. We take the parameter 7 to be real, which corresponds 
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FIG. 3: The line shapes in a short-distance channel (solid lines) and in the D^D^tt^ channel (dotted 
lines). The line shapes are shown for five values of 7: —48, —34, 0, 34, and 48 MeV. The line shapes 
in the short-distance channel and in the D^D^ir^ channel are separately normalized so that the 
curve for 7 = has a maximum value of 1. Kl E = 0.2 MeV, the order of both the solid curves 
and the dashed curves from top to bottom is 7 = 0, 34, —34, 48, and —48 MeV. 



to the assumption that decay modes other than D°£)°7r° and D°D°7 give a small contribution 

to the total width of X(3872). In Fig. 3, the solid lines are the line shapes for a short-distance 
decay mode, which is given by Eq. (26). We neglect the energy-dependence from the factor 
T'-'^E), which is very small for J/ipT^^n^ . The dashed lines are the line shapes for D'^D^ir^, 
which is given by Eq. (27). We show the line shapes for five values of 7: 48, 34, 0, —34, 
and —48 MeV. For 7 = 34, 0, and 48 MeV, the peaks of the resonance are at = —0.6, 
0, and —1.2 MeV, respectively. They correspond to the central value of the measurement 
in Eq. (1) and deviations by ilcr from the central value. In Fig. 3, the line shapes for the 
short-distance channel and for the D'^D^n^ channel are separately normahzed so that the 
maximum value is 1 for 7 = 0. This figure illustrates that as 7 decreases toward from 
above, the area under the short-distance line shape decreases and it becomes very small 
for negative values of 7. As 7 decreases toward from above, the area under the D^D^tt'^ 
line shape decreases less rapidly, but it is also small for negative values of 7. At the values 
7 = 48, 34, 0, —34, and —48 MeV, the differences between the positions of the peaks in the 
D^D^-K^ channel and the short-distance channel are 0.00024, 0.00052, 0.019, 0.75, and 1.75 
MeV, respectively. Thus the 4.7^}'° MeV difference between the masses in Eqs. (2) and (1) 
appears to be only compatible with negative values of 7. 

In an actual measurement, the line shapes shown in Fig. 3 would be smeared by the effects 
of experimental resolution. In the Belle discovery paper, the X(3872) signal in J /ipT:^T:~ was 
fit by a Gaussian with width 2.5 MeV, which is compatible with the experimental resolution 
[1]. For the D^D^tt^ enhancement near threshold, the signal was fit by a Gaussian with 
width 2.24 MeV, which is again compatible with the experimental resolution [5]. In Fig. 4, 
we illustrate the effects of experimental resolution by smearing the line shapes in Fig. 3 
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E (MeV) 

FIG. 4: Same as in Fig. 3, except that the hne shapes have been smeared by a Gaussian function 
with width 2.5 MeV to mimic the effects of experimental resolution. The smeared line shapes in 
the short-distance channel and in the D^D^tt^ channel are separately normalized so that the curve 
for 7 = has a maximum value of 1. At -E = MeV, the order of both the solid curves and the 
dashed curves from top to bottom is 7 = 48, 34, 0, —34, and —48 MeV. 

by a Gaussian with width 2.5 MeV. At the values 7 =_48, 34, 0, -34, and -48 MeV, the 
differences between the positions of the peaks in the D^D^n^ channel and the short-distance 

channel are 0.45, 0.82, 2.0, 2.4, and 3.2 MeV, respectively. Thus, after taking into account 
the effects of experimental resolution, the 4.7lj^'g MeV difference between the masses in 
Eqs. (2) and (1) is not incompatible with a positive value of 7. 



V. ESTIMATE OF THE SHORT-DISTANCE FACTOR 

In this section, we give an order-of-magnitude estimate of the short-distance factor 
that appears in Eqs. (25), (26), and (27). We use the measurement by the Babar Collabo- 
ration of the branching fraction for the decay of into K~^D*^D^ [26]: 

Br[S+ ^ X+D*°L>°] = (4.7 ± 1.0) x 10"^ (29) 

The upper bound on the branching fraction for K^D^D*^ is 3.7 X 10-^ at the 90% confidence 
level. The partial width for the decay in Eq. (29) can be written 

r[s+ ^ x+l>*ol>°] = ^ (|A^r>*f+^*°'^°, (30) 

1Mb 

where c^^^^*^^^ = 194.5 MeV^ is the integral of the three-body phase space and 3(|A^p) is 
the square of the matrix element averaged over the relativistic three-body phase space and 
summed over the D*° spins. Using the measured hfetime, t[S+] = 1.64 x 10~^^ s, and the 
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central value of the branching fraction in Eq. (29), we find that the average value of the 
square of the matrix element in Eq. (30) is = 3.4 x 10~^^. 

The differential rate for to decay into K+ + D*^D^ with D*^D^ near its threshold can 
be written 



[B. ^ K- + D-'D"] = '""'^"^l'':^ ''"'^•'^ \M?V^, (31) 



dV 

— ... ..... J - 647r3M3 



where E is the rest energy of the D*^D^ relative to its threshold and is the square 

of the matrix element averaged over angles in the D*^D^ rest frame and averaged over the 
D*^ spins. The dramatic dependence of \}A\^ on E comes from the resonance factor \ f[E)\'^ 
in Eq. (27). By comparing Eq. (31) with Eq. (27), we can obtain an expression for the 
short-distance factor r^+: 

_ 3AV^(M^,M.o + Mo,m^) ^ , 
Y28^^ eV^\ ■ 

Because of the threshold resonance, the square of the matrix clement |A^p for — > 
+ D*^ should be significantly larger near the D*^D^ threshold than its average 
over the entire three-body phase. We will assume that this threshold resonance is the only 
effect that gives a significant enhancement of in this corner of the phase space. Thus 
as an order of magnitude estimate of |A^P, we will take 

i^i^-^i^i^^fttW- '''' 

where is the average value of the square of the matrix element defined by Eq. (30) 

and A is the momentum scale below which the resonant behavior sets in. The natural scale 
for A is (2M,iii/)V2 = 125 MeV or = 135 MeV. The estimate in Eq. (33) is hkely to be 
an underestimate because is enhanced by various resonant contributions, including 

Dsi{2536) + -D°. Using the value of that corresponds to the central value of the 

branching fraction in Eq. (29), our estimate of the short-distance factor is 

A2 

rf+ « (3.8 X 10^^^ MeV) — . (34) 



VI. ANALYSIS OF REF. [25] 



In Ref. [25], the authors analyzed the data from the Belle and Babar collaborations on 
B ^ K + X(3872) using a generahzed Flatte parameterization for the {D*D)^ scattering 
amplitude: 

/hkME) =(^-^[E-Ef + iT^+,-j/^{E)/2 + iT,+,-,oj,^{E)/2\ + k{E) + k^{E)^ , 

(35) 

where k{E) = (-2M,oo^ + ^e)^/^ and Ki{E) = {-2M^ii{E - + lef/^. The functions 
T^+^-jl^{E) and Tt^+„-^oji^{E) are determined by the p and uj resonance parameters up to 
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normalization factors fp and f^^. The resonant term in the differential branching fraction for 
the inclusive decay of into was expressed in the form 

dBv 2 

— [B+ -^K+ + resonant] = -{B/g) lmf{E), (36) 
dE TT 

The adjustable parameters in the model of Ref. [25] are g, Ef, fp, f^, and B. The imaginary 
part of the scattering amplitude in Eq. (35) is 

Im/HKKN(^) = |/(^)r [T,+,-j/^{E)/g + V^+^-,,jj^{E)/g - lmK{E) - Im«;i(E)] . 

(37) 

The first two terms in the square brackets are identified with the contributions from 
J/-^ TT+TT" and J/V' tt"'"tt~tt'', respectively. The third term —\m.K{E^ is nonzero only for 
£■ > and was identified with the contribution from D^D'^n^. The last term — Imfi;i(£') 
is nonzero only ior E > u and could be identified with the contributions from D*^D^ and 

In the analysis of Ref. [25], the ratio fuj/fp was determined from the branching ratio 
of X(3872) into J/ip n^n'n^ and J/ip7i~^n~ measured by the Belle collaboration [5]. The 
parameters g, Ef, fp, and B were then used to fit the J/z/'tt+tt" energy distributions in the 
decay — > K'^ + J/ipn^ii^ measured by the Belle and Babar Collaborations [1, 27] and 
the D^D^n^ energy distribution in the decay B~^ — > + D^D^ti^ measured by the Belle 
collaboration [8]. The authors of Ref. [25] found that their fits exhibited scaling behavior. 
The fits were insensitive to rescaling all the parameters g, Ef, fp, f^, and by a common 
factor. From the expressions in Eqs. (35), (36), and (37), one can see that the scaling 
behavior simply implies that the term —{2/g)E in the denominator of Eq. (35) can be 
neglected. The authors interpreted the parameter gf as a DD* coupling constant. However 
a more appropriate interpretation of g is that it determines the effective range for S-wave 
scattering in the {D*D)^ channel. If we neglect the small imaginary terms proportional 
to r^+^-j/^(£J) and r^+^-7rOj/^(£') in the denominator in Eq. (35), the inverse scattering 
length and the effective range are 

7 ^ -^{Ef/g) - ^/2M^, (38a) 
rs ^ -2/(M,oo^?)-l/\/2AW. (38b) 

We have simplified the expression for in Eq. (38b) by setting M*ii/M^,oo = 1- The 
inverse scattering lengths 7 calculated from Eq. (38a) using the parameters of the fits in 
Ref. [25] range from —51.7 MeV to —66.3 MeV. These values are close to the exact results 
for the complex inverse scattering length 7 for the fits in Ref. [25]. Their real parts ranged 
from —48.9 MeV to —63.5 MeV and their imaginary parts ranged from 3.3 MeV to 6.8 
MeV. Because of the scaling behavior, the scattering amplitude near the D*^D^ threshold 
is insensitive to the effective range r^, so its numerical value cannot be determined from the 
fits of Ref. [25]. The authors could have obtained essentially the same fits by taking the 
range to be 0. They could have replaced the scattering amphtude in Eq. (35) by 

/hkkn(^) ~ [-Re 7 - iT^+^-ji^{E)/g - ir^+^-^oj/^{E)/g + k{E)]~^ , (39) 
with Re 7, fp/g, f^/g, and B/g as the adjustable parameters. 
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The authors of Ref. [25] found that they could obtain acceptable fits to the Belle and 
Babar data only if the J/ip-K^Ti' energy distribution was peaked almost exactly at the 
]j*ojjo threshold and the D^D^n^ energy distribution was peaked 2 to 3 MeV above the 
£)*o£)0 threshold. This requires a negative value for Re 7, which would correspond to a 
virtual state. They concluded that if the D^D^t:^ threshold enhancement is associated with 
the X(3872), then X must be a virtual state of charm mesons. 

A crucial flaw in the analysis of Ref. [25] is that they assumed that the D^D^tt^ contri- 
bution was proportional to —Im /«(£'), where k{E) = {—2M^ooE — ieY^^. This expression 
vanishes if £■ < 0. If the X(3872) is a bound state with Re7 > 0, most of the support of 
its line shape is in the region < 0. Thus the assumption of Ref. [25] essentially forbids 
the bound state from decaying into D^D^n^. This is in direct contradiction to one of the 
universal features of an S-wave resonance near threshold. As the scattering length a in- 
creases to -|-oo, the binding energy approaches l/(M^,oo o^) and the mean separation of the 
constituents approaches a/2. In this limit, the decay of the resonance should be dominated 
by decays of its constituents. Since D^n^ is the largest decay mode of D*^, D^D^n^ should 
be a dominant decay mode of the X(3872) if it is a bound state. This flaw in the analysis of 
Ref. [25] could be removed by using the expression for k{E) in Eq. (20) to take into account 
the effects of the decays of the constituents of the X(3872) resonance. 

By comparing Eqs. (25) and (36), we can express the short-distance factor in terms 
of the parameters of Ref. [25] : 

rS^-iB/g)T[B^]. (40) 

TT 

The best fits in Ref. [25] give values of rf+ that range from 3.8 x lO^^^ y^^y to 4.9 x 10"^^ 
MeV. This is one order of magnitude larger than the estimate for the short-distance factor 
in Eq. (34) with A = rrij^. This requires the square of the matrix element for — > 

+ D*^D^ near the D*^D^ threshold to be larger than the average value defined 
by Eq. (30) not only by the threshold resonance factor m^|/(£')p but also by an additional 
order of magnitude. 

VII. FITS TO THE ENERGY DISTRIBUTIONS 

The line shapes of A(3872) in the J/ipn^n' and D^D^n^ decay channels have been 
measured by the Belle and Babar Collaborations for the production processes B ^ K + X. 
In this section, we fit the Belle measurements for the production process S+ — > K'^ -\- X to 
a theoretical model for the line shapes that takes into account the width of the constituent 
D*^ as well as inelastic scattering channels for the charm mesons. 

A. Experimental data 

The Belle and Babar Collaborations have both measured the energy distribution of 
J/V'TT+TT- near the A:(3872) resonance for the decay B+ ^ + J/iI^-k+'k- [1, 4]. The 
Babar data has larger error bars, so we will consider only the Belle data. The Belle data 
on the J/ipTT^Ti' energy distribution is given in Fig. 2b of Ref. [1]. The figure shows the 
number of events per 5 MeV bin as a function of Mj/^^+„- from 3820 MeV to 3920 MeV. If 
we use the CLEO measurement of the mass and the PDG value for the D*^ — D° mass 
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E N AN 



—19.3 


—2.10 


2.78 


—14.3 


— 1.10 


2.89 


— 9.3 


— 1.21 


3.04 


-4.3 


9.60 


4.83 


0.7 


24.56 


6.46 


5.7 


-1.47 


2.99 


10.7 


-1.57 


2.99 


15.7 


-0.58 


3.25 


last 7 bins 


28.23 


10.54 



TABLE I: Belle data on the J/V'tt+tt" energy distribution: numbers of events N and their uncer- 
tainties AN in 5 MeV bins centered at the energies E. The numbers of events A'^ were obtained 
from Fig. 2b of Ref. [1] by subtracting the linear experimental background. We used only the last 
7 data points in our analysis. 

difference, this corresponds to E extending from —51.8 MeV to +48.2 MeV. In Ref. [25], 
Hanhart et al. used only the 8 bins extending from —21.8 MeV to +18.2 MeV. They sub- 
tracted the linear experimental background shown in Fig. 2b of Ref. [1] to get the 8 data 
points given in Table I. Our values for the energies at the centers of the bins differ from 
those in Ref. [25] by —0.2 MeV, because they used the PDG values for the masses Mq and 
M*o instead of using the CLEO value for Mq and the PDG value for M*o — Mq. In our 
analysis, we have chosen to omit the first data point in Table I so that the data points are 
more symmetric about E = 0. Wc will use a theoretical model for the energy distribution 
that is essentially at = —19.3 MeV, so including this point would simply increase the 
by a constant 0.57. The 7 data points used in our analysis are plotted in Fig. 5. Note 
that there are only two bins in which the data differs from by significantly more than one 
error bar. 



E 


A 


AN 


-5.015 


0.42 


1.03 


-0.765 


0.90 


1.38 


3.485 


11.58 


4.13 


7.735 


1.35 


2.82 


11.985 


1.50 


3.30 


16.235 


-0.89 


3.09 


all 6 bins 


14.86 


6.96 



TABLE II: Belle data on the L>"L»" TT energy distribution: numbers of events N and their uncer- 
tainties AN in 4.25 MeV bins centered at the energies E. The numbers of events N were obtained 
from Ref. [28] by subtracting the total experimental background. Only the last 5 data points were 
used in the analysis of Ref. [25] . 

The Belle Collaboration has measured the energy distribution of J/-?/' 7r"'"7r~ near the 
X(3872) resonance for the decay B ^ K + D^D^ti^ [8]. The D^D^tt^ energy distribution 
is shown in Fig. 2a of Ref. [8]. The figure shows the events per 4.25 MeV bin as a function 
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FIG. 5: Number of events per 5 MeV bin for + J/ipir~^Tr~ as a function of the total 

energy E of J/'ip7r~^-K~ relative to the D*^D^ threshold. The data points, which are given in 
Table I, were obtained by subtracting the linear experimental background from the data of the 
Belle collaboration in Ref. [1]. The theoretical curves are the differential number distributions 
dN/dE in Eq. (43) multiplied by the 5 MeV bin width. The three curves correspond to the global 
minimum of (dashed curve), the local minimum of (solid curve), and the fit ^Belle of R-^f. [25] 
(dotted curve). The inset shows the peaks of the distributions. 

of MjyoQo^o - 2Mdo - M^o from MeV to 76.5 MeV. This corresponds to E extending 
from -7.14 MeV to +69.36 MeV. In Fig. 2a of Ref. [8], the data for B+ ^ K+ + X and 
+ X are combined in the same plot. The energy distributions for decay and 
B^ separately were presented at the ICHEP 2006 conference [28]. In Ref. [25], Hanhart et 
al. used only the data for B^ — > + X m b bins extending from —2.89 MeV to +18.36 
MeV. They subtracted the combinatorial background to obtained the data points for those 5 
bins. To account for the remaining experimental background, which is an increasing function 
of -E, they added a background term to the theoretical expression for dV/dE and determined 
its coefficient by fitting to the data. That background term is weakly constrained by the 5 
data points. We have therefore chosen to subtract the total experimental background instead 
of only the combinatorial background. The resulting data points are given in Table II. Our 
values for the energies at the centers of the bins differ from those in Ref. [25] by —0.34, 
because they used the PDG values for the masses Mq and M^,q to determine the D*^D^ 
threshold instead of using the CLEO value for Mq and the PDG value for M^q — Mq. The 
first data point in Table II was omitted in the analysis of Ref. [25] , because it would have 
given a constant contribution to of 0.17. In our analysis, we have chosen to include this 
data point even though its effect on our analysis is negligible. The 6 data points used in our 
analysis are plotted in Fig. 6. Note that there is only one bin in which the data differs from 
by significantly more than one error bar. 

The analysis of Ref. [25] also used the Belle measurement of the branching ratio for decays 
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FIG. 6: Number of events per 4.25 MeV bin for B+ ^ K+ + £)°Z)°7r° as a function of the 
total energy E of D^D^tt^ relative to the D*^D^ threshold. The data points, which arc given in 
Table II, were obtained by subtracting the total experimental background from the data of the 
Belle Collaboration in Ref. [8]. The theoretical curves are the differential number distributions 
dN/dE in Eq. (44) multiplied by the 4.25 MeV bin width. The three curves correspond to the 
global minimum of (dashed curve), the local minimum of (solid curve), and the fit ^Beiie of 
Ref. [25] (dotted curve). 



of X (3872) into J/V^ 7r+7r-7r° and J/V^tt+tt" [5]: 

Bt[X J/V>7r+7r-7r°] 



1.0 ±0.4 ±0.3. 



(41) 



The signal region for J/ipir^n 7r° included only energies E within 16.5 MeV of 3872 MeV. 
Using the CLEO value for Mq and the PDG value for M*o — Mq, this corresponds to 
-16.3 MeV <E < +16.7 MeV. 



B. Theoretical model 



We summarize our model for the X(3872) lines shapes in the channels J/ip 7r~^7:~ , 
J/'07r+7r~7r°, and D^D^n^ channels. The J/'^tt+tt" and J/V' 7r+7r~7r° energy distributions 
are given by Eq. (26), while the D^D^t^^ energy distribution is given in Eq. (27). The scatter- 
ing amplitude f{E) is given by Eq. (16) with Im7 replaced by V^l^^^^' {E)+r'^/'t"'^''~''° (E): 



f{E) = [-Re 7 - iV-^/-^^^^ (E) - iT 



(E) + k{E) 



(42) 



The energy dependence of the functions T-^/^^^^ {E) and T'^/^^^^ ^°{E) are shown as solid 
hues in Fig. 2. Their normalizations are determined by tJ'^'' = r'^/i'^+^- (^Q) and r^^"" = 
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^"(0), which we treat as adjustable parameters. The function k{E) in Eq. (42) 
is given in Eq. (20). The functions r*o(-E) and Brooo(-E') are given in Eqs. (11) and (12), 
respectively. Thus our model has 4 adjustable real parameters: Re 7, r^+, Fq^'^, and Tq^^. 

To translate the differential rates dV/dE into numbers of events in the Belle experiment, 
we follow the prescription used in Ref. [25]. The differential number of J/'^tt+tt" events at 
the energy E is 

^rjM,+,-l ^ N^SJ r[B-]/h 

^^[^/^TT TT J B^^B+^j^+^x]Br[X^J/iJ7r+7r-]dE^'' " +J/^7r7rJ. 

(43) 

For the number of observed events, we use the central value in Ref. [1]: -^.^strved ~ 35.7. 
For the product of branching fractions in the denominator, we use the central value in 
Ref. [1], which is 1.3 x 10~^. The differential number of D^D^tt^ events at the energy E is 

— fD°L)°7r°l = ^ibsorvSi r[B-]/h ^ ^ ^+ D°L)°7r°l . (44) 

dE^ ^ Bt[B+ ^ K+ + D0D^7r0]dE^ ^ ^ ' 

For the number of observed events, we use the central value in Ref. [8]: iV^^L^vcd = 17.4. 
For the branching fraction in the denominator, we use the central value in Ref. [8] , which is 
1.02 X 10^1 

Following Ref. [25] , we take into account the Belle measurement of the branching ratio in 
Eq. (41) through a constraint on the parameters of our theoretical model. We demand that 

1.0. (45) 

This constraint determines the ratio Fq^'^/Fq^'^ as a function of the parameters Re 7 and 
pV'27r j£ combine the errors in Eq. (41) in quadrature, the error bar on the right side of 
Eq. (45) is ±0.5. We ignore this error bar and constrain the ratio in Eq. (45) to be 1.0 so 
that our results can be compared more directly with those in Ref. [25].^ 



r[B+- 


> K+ + J/'^|;7^+7^- 


7r°;|£;| < 16.5 MeV] 


r[B+ - 


^ K+ + J/i/jn+n 


-;\E\< 16.5 MeV] 



C. Fitting procedure 

One of the most important issues is whether the data on the X(3872) is compatible with 
it being a bound state (corresponding to Re 7 > 0) or whether it must be a virtual state 
(corresponding to Re 7 < 0) as advocated in Ref. [25] . We therefore analyze the Belle data 
by fixing the parameter Re 7 and minimizing the with respect to the other 3 adjustable 
parameters Fg+ , Fq^'^, and Fq^'^ subject to the constraint in Eq. (45). The is the sum of 
13 terms corresponding to the last 7 data points in Table I and the 6 data points in Table II. 
The constraint in Eq. (45) determines the ratio F^^'^/F^^'^ for fixed Re 7. In Fig. 7, we show 
the minimum value of with respect to variations of the two remaining parameters F^+ 
and Fq^'^ as a function of Re 7. The global minimum is = 8-3 at Re 7 = 57.8 MeV. There 
is also a local minimum at Re 7 = 17.3 MeV with = 10.0. Since Re7 > for both the 



^ In Ref. [25], the partial widths on the right side of Eq. (45) were integrated over the larger region \E\ < 20 
MeV, but this difference has a negligible effect on the analysis. 
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ReY(MeV) 

FIG. 7: The minimum for the 13 data points in Figs. 5 and 6 as a function of Re 7. The has 
been minimized with respect to the parameters , Fq^^, and Fq^'^ subject to the constraint in 
Eq. (45). The dot is the for the fit ^Beiie of Ref. [25]. 



fit 


Re7 


Ff + X 


-pV'27r 
^ 


^ 




global minimum 


+57.8 


4.1 


0.44 


0.50 


8.3 


local minimum 


+17.3 


8.3 


2.3 


2.4 


10.0 


^Belle of Ref. [25] 


-48.9 


38 


0.93 


0.71 


19.1 



TABLE III: Parameters (in units of MeV) and values of x^ for three fits to the subtracted Belle 
data: the global minimum of x^, the local minimum of x^, and the fit ^Beiie of Ref. [25]. 

global minimum and the local minimum, these fits correspond to a bound state. For the 
sake of comparison, the AbcUc fit of Ref. [25] gives x^ = 19.1. The real part of the inverse 
scattering length for this fit is Re 7 = —48.9 MeV. Since this is negative, this fit corresponds 
to a virtual state. In Fig. 7, the value of for the ^Beiie fit is shown as a dot that lies just 
above the fine. 

In Table III, we list the parameters and the values of x^ for the global minimum, the 
local minimum, and the fit A-qcWc of Ref. [25]. For the fit AbeIIo R-e7 is the real part of 
the inverse scattering length, is given in Eq. (40), and Vf^ and Fq^'^ are the values of 
TT,+T,-jj^{E)/g and r^+^-^oj/^{E)/g a.t E = 0. Our order-of-magnitude estimate for F^+ is 
given in Eq. (34) with A = m„. The values of for the global minimum and the local 
minimum of x^ are more consistent with this estimate than the fit AbcUc of Ref. [25]. 

The line shapes of X(3872) in the J/ip 7r~^7r~ and D^D^n^ decay channels for various 
fits are shown together with the Belle data in Figs. 5 and 6, respectively. The line shapes 
corresponding to the local minimum, the global minimum, and the fit ^Beiie of Ref. [25] 
are shown as solid, dashed, and dotted lines, respectively. The data points in these figures 
should be compared to the average values of the curves over the appropriate bins centered 
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on the data points. The global minimum of is somewhat pathological in that for both 
J/%l)TT^TT~ and D^D^Ti^ the integral of the line shape over a bin is largest not in the bin 
with the highest data point, but in the next lower bin. For the local minimum of x^, the 
integrated line shape gives a good fit to the data point in the highest bin for J/ipTr~^Tr~ but 
it is smaller than the data point in the highest bin for D^D^n^ by more than two standard 
deviations. 

VIII. DISCUSSION 

We have developed an approximation to the line shapes of the X(3872) that takes into 
account the width of its constituents D*^ or D*^ as well as inelastic scattering channels for 
D*^D^ and D°L>*°. The best combined fit to the Belle data on the energy distributions 
for the X resonance in the J/iI)ti^ti~ and D^D^ti^ channels corresponds to a bound state 
with Re 7 > 0, although a virtual state with Re 7 < is not excluded. Our results are 
in contradiction to the conclusions of Ref. [25]. By fitting essentially the same data, they 
concluded that the X(3872) must be a virtual state. The flaw in their analysis was that 
they did not allow for decays of the constituent D*° or D*^ in the case where the X is a 
bound state. 

Our theoretical model for the line shapes is based on the assumption that the inverse 
scattering length in the {D*D)^_^ channel is small compared to all other relevant momentum 
scales. The line shapes are determined by the scattering amplitude f{E) for charm mesons 
m the {D*D)l channel given in Eq. (42). This amplitude contains no information about 
some of the other nearby thresholds, including the D*^D~ threshold at £^ = +8.1 MeV 
and the £)°L>°7r° threshold at £^ = —7.1 MeV. Thus the hne shapes can be expected to be 
accurate only when \E\ is much less than 7 MeV. However the maxima of the line shapes 
are well within the region of validity. For the local minimum of x^, the peaks in the line 
shapes are near E = —0.17 MeV for J/ijj7r^7r~ and near E = —0.14 MeV for D^D^tt^. For 
the global minimum of x^, the peaks in the line shapes are near E — —1.7 MeV for both 
J/'^TT+TT" and D^D^n^. The effects of the other thresholds would be dramatic only when 
is large enough that the line shapes are already small. They would be unlikely to change 
our qualitative conclusion that the best fit to the data corresponds to a bound state. 

The range of validity of our model for the line shapes could be extended by taking into 
account exphcitly scattering channels for the charged charm mesons D*^D~ and D^D*~. 
A first step in this direction has been taken by Voloshin [29]. He showed that the isospin 
symmetry breaking pattern of QCD provides interesting constraints on the line shapes. In 
particular, he predicted a zero in the line shape of X(3872) in the J/ipTT~^7:~7i^ channel 
between the D*^D^ and D*^D~ thresholds. We have carried out a more rigorous two- 
channel analysis [32]. Given a plausible dynamical assumption, we find that the line shape 
of X in the J/ipTT^n'n^ channel has a zero between the D*^D^ and D*^D~ thresholds in 
the decays + X but not in the decays + X. We also find that the 

line shape of X in the J/'^tt+tt" channel can have a zero below the D*^D^ threshold in the 
decays B'^^K'^ + X. 

The accuracy of our predictions for the line shapes could be further improved by taking 
into account pions explicitly. The system consisting of D*^D^, D^D*^ ^ and D^D^n^ states 
with energies near the D^D*^ threshold can be described by a nonrelativistic effective field 
theory. The simplest such theory has S-wave scattering in the {D*D)^_^ channel and 
couphngs that allow the decay D*° — > D^n^. Fleming, Kusunoki, Mehen, and van Kolck 
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developed power-counting nilcs for this effective field theory and showed that the pion 
couplings can be treated perturbatively [30]. They used the effective field theory to calculate 
the decay rate for X(3872) — > D°Z)°7r° to next-to-leading order in the pion coupling. 

In applying this effective field theory to the line shapes of the X(3872), one comphcation 
that will be encountered is infrared singularities at the D*^D^ threshold that are related 
to the decay D*^ D^n^. This problem has been analyzed in a simpler model with spin- 
particles and momentum- independent interactions [31]. The problem was solved by a 
resummation of perturbation theory that takes into account the perturbative shift of the 
]jo£)*o threshold into the complex energy plane because of the nonzero width of the D*^. 

In summary, the establishment of the quantum numbers of the X{3S72) as 1"^"*" and the 
measurement of its mass imply that it is either a charm meson molecule or a charm meson 
virtual state. The existing data favor a charm meson molecule, but a virtual state is not 
excluded. To decide conclusively between these two possibilities will require more extensive 
data on the line shapes of the X{3872) in various decay channels and for various production 
processes. 
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